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The article analyzes similarity of closed polygonal curves with respect to the Frechet 
metric, which is stronger than the well-known Hausdorff metric and therefore is more 
appropriate in some applications. An algorithm is described that determines whether 
the Prechet distance between two closed polygonal curves with m and n vertices is less 
than a given number e. The algorithm takes 0{mn) time whereas the previously known 
algorithms take 0{mn\og{rnn)) time. 

Keywords: computational geometry, Frechet distance, computational complexity. 


1. Introduction 

The Frechet metric is used for cyclic process analysis and image processing [5]. It 
is stronger than the well-known Hausdorff metric 0, i and therefore is more 
appropriate in some applications [T] . The Frechet metric for closed polygonal curves 
has been studied in a paper [T] by Alt and Godau. They propose an algorithm that 
determines whether the distance between two closed polygonal curves with m and 
n vertices is greater than a given number e. The complexity of the algorithm is 
0{mnlog{mn)) on a random access machine that performs arithmetical operations 
and computes square roots in constant time. Our paper shows that the compu- 
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tational complexity of the problem is less than 0{mnlog(rnn)) and provides an 
algorithm that takes 0(mn) time to solve the problem. The exact formulation of 
the problem is given in Section [2l Sections |3] describes the concepts of the original 
paper [T], which with slight modifications serve as the basis for our paper. The 
difference between the proposed and known approaches is specified at the end of 
Section [3l Sections |4][6] describe the proposed approach. 


2. Problem definition. 

Let be a linear space with the Euclidean distance d : x K.^' —)• K. 

Definition 1. A closed m-gonal curve X is a pair {x,fx) where x is a sequence 
{xq, Xi, - ■ ■ , Xm = ;co), Xi S and fx is a function [ 0 , m] ^ such that 
fxif + a) = (1 — a)xi + axi+i for i £ {0,1,..., m — 1} and 0 < a < 1. 

Definition 2. A cyclic shift of an interval [0, m] by a value t £ [0, m] is a function 
s : [0,to] [0,m.] that depends on a parameter t such that s{t\T) = t + r for 

t + T < m and s(t; r) = t + t — m for t + t > m. 

Eor any number m let Wm be the set of all monotonically non-decreasing con¬ 
tinuous functions [0,1] ^ [0,m] such that w(0) = 0, w(l) = m. 

Definition 3. A function ip : [0, 1] —>■ is called a monotone reparametrization of 
a closed m-gonal curve X = {x, fx) if a function w £ Wm and a number r £ [0, m] 
exist such that ip{t) = fx{s{w{t);T)) for all t £ [0,1]. 

Eor given closed polygonal curves X and Y denote ‘I’x and sets of their 
reparametrizations. 

Definition 4. The Frechet distance between closed polygonal curves X and Y is 
6{X,Y)= min min max (t),(^Y(t)). 

0<i<l 

The problem consists in developing an algorithm that determines whether 
6(X,Y) < e for given closed polygonal curves X and Y and a number e. 


3. The free space diagram and pointers. 


The problem’s analysis is based on the concept of a free space diagram introduced 
by Alt and Godau [T] in the following way. For two numbers m and n let us define 
a rectangle D = [0,m] x [0,n] with points (u,u) £ D. For two closed polygonal 
curves X and Y with m and n vertices and a number e a subset Dg = {(u,u) £ 


D 


difxiu), friv)) < e} is defined. Let us also define a rectangle D = DU {(m ■ 


m, v) {u, v) £ D} with its subset Dg = DgU{{u-\-m, v) {u, v) G Dg} called a free 
space. Denote T, B, L and R the top, bottom, left and right sides of the rectangle 
D. Denote Dij a subset [i — l,f] x [j — l,j] and call it a cell. 
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Definition 5. A monotone non-decreasing path (or simply, a monotone path) is 
a connected subset 7 C such that [u — u'){v — ^ 0 fof any two points 

{u,v) e 7, {u',v') G 7. 

Note that this definition allows a monotone path to contain vertical segments. 
That is why a condition {u — u'){v — v') > 0 is used instead of standard form 
(u — v')/(u — u') > 0 of the definition of a non-decreasing function v = f{u). 

Definition 6. Two points {u,v) G D and (u',v') G D are mutually reachable if and 
only if a monotone path 7 exists such that {u,v) G 7, (u',v') G 7. 

Definition 7. A point {u, v) G is reachable from the bottom if it is reachable 
from at least one point from B; a point (u, v) G Dg is reachable from the top if it is 
reachable from at least one point of T. 

Denote 54 . C Dg a set of points reachable from the bottom and C a set of 
points reachable from the top. 

Let us define two pointer functions —>■ [0, 2m] and ?' 4 , : —>■ [0, 2m]. 

For (u,v) G the pointer r'^(u,v) is the maximum value u* such that {u,v) is 
reachable from (u*,n) G T. For {u, v) € gi\ B the pointer ri(u, v) is the maximum 
value u* such that (it, v) is reachable from (it*, 0) G B. For (it, 0) G gfCiB the pointer 
r 4 ,(it, 0) equals it. Figure [2] illustrates the introduced concepts and designations. 



Fig. 1: Two closed polygonal curves with an interval of length e 


We rely on the following lemma proved in the paper [T] by Alt and Godau (see 
lemma 9 m)- 

Lemma 1. The distance between closed polygonal curves X and Y is not greater 
than e if and only if there exists a number u G [0, m], such that the points (u-fm, n) 
and (it, 0) are mutually reachable. 

The following lemma is similar to Lemma 10 [T] as well as its proof. 
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Fig. 2: The free space diagrani. The light area is the set De- 
r"^{u 2 ,V 2 ) is the rightmost reachable point on T from {u 2 ,V 2 )- 
fliuijVi) is the rightmost reachable point on B from (ui,ni). 


Lemma 2. Two points {ut,n) G T and (ub,0) G B are mutually reachable if and 
only if 

{ut,n) € gi, {uh,Q)€g^, m* < r'l'(ub, 0), Ub<r^{ut,n). 

Proof. Obviously, if {ut,n) G T and {ub,0) G B are mutually reachable then 
(ut, n) G gi, {ub, 0) G g'^ and Ut < 0), Ub < ri{ut, n). 

The reverse implication is also valid, which is illustrated by Figure |3l Let ■yt C 
be a monotone path from {ut,n) to {ri{ut,n),0) and yb C be a monotone path 
from {ub, 0) to {r'^{ub, 0), n). Both paths are connected subsets that due to conditions 
Ut < r^{ub,0), Ub < r^_{ut,n) intersect in at least one point {uq,vq). Let us consider 
a path 7 that consists of a segment of yb from {ub,0) to {uq,vo) and a segment of 
7 * from (uojVq) to {ut,n). The path y is monotone, it is contained inside and 
connects (ub,0) and {ut,n). □ 

According to Lemmas [Hand [2] testing condition S{X, T) < e is reduced to finding 
value u that fulfills 

(u,0 )Gg^, {u + m,n)&gi, M + m<r^(u,0 ), u < r^iu + m^n). (1) 

The pointer functions r'^ and are similar to the pointers defined by Alt and 
Godau. However, the pointer function takes values from T and function takes 
values from B, while Alt and Godau consider pointers with values from T L) R and 
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with values from L\J B. The Alt’s and Godau’s pointers allow to use a divide and 
conquer type algorithm that merges either vertically or horizontally two diagrams 
with known pointers and obtains a bigger diagram with pointers for the new dia¬ 
gram. The pointers can be computed in 0(1) time for diagrams containing only one 
cell. By sequentially merging smaller diagrams into bigger ones, the pointers for the 
whole diagram can be obtained in 0(rnnlog{mn)) time. 

In this paper we rely on a recurrent relation between pointer values and 
on cell borders. It is trivial to compute on T and on B. Our algorithm does 
not use the divide and conquer approach but proceeds cell by cell. It propagates 
pointers on each cell’s borders using the recurrent relation and eventually obtains 
the values on B and on T in 0{mn) time. 

4. Formal properties of pointer functions. 

For each 1 < i < 2m and 1 < j < n denote and Rij the top and right borders 
of a square cell Dij = [i — l,i] x [j — l,j] and extend these denotations so that Tio 
is a bottom border of Du and Rqj is a left border of Dij. Let us designate 

Bij = Lij = 

TRij = Rj U Rij , LBij = Lij U B,j . 

In order to test 0 the sets gi and functions have to be expressed with 

a finite data structure. It is significant that intersection n Dij is convex for each 
pair (i,j) [1]. It is not difficult to prove that for each pair {i,j) the intersections 
gi n Tij, g^ n Rij are also convex, so that each of these intersections is a single 
interval (line segment) on the border of Dij. Simple recurrent relations hold for 
these intervals, so that intervals g^ H and g^j n Rij can be computed based on 

gi n Bij and g_i fl Lij as well as fl Bij and gb fl Lij can be computed based on 
g'^ n Tij and g"^ fl Rij. The Algorithm 1 described in [T] for slightly different purposes 
is an example of this computation. Each step of the recursion takes constant time. 
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consequently, computing intersections fl T^-, fl Rij , fl Bij and g^ fl for 
all {i,j) takes 0(mn) time. Therefore, from now on it is assumed that these line 
segments are available. 

As for the functions r'^ and they require more detailed analysis. Let us de- 
hne partial ordering ^ on the set D such that (Mi,r;i) ^ (u 2 ,V 2 ) if and only if 
ui < U 2 and vi > V 2 - On each subset TRij and LBij relation ^ is a complete 
ordering. Therefore, for each closed subset b C TRij and for each closed subset 
b C LBij a symbol extrfe will be used as a designation of such point {u*,v*) G b 

that {u',v') =4 {u*,v*) for each {u',v') G b. Both and r"! are monotone functions 

of their argument in a sense of the following lemma. 

Lemma 3. 

Let {u 2 ,V 2 ) G gi- If {ui,vi) 4 (u 2 ,U 2 ) then ri{ui,vi) < ri{u 2 ,V 2 ). 

Let {u2,V2) G g^. If (ui,vi) ^ {u2,V2) then rl{ui,vi) < rI{u2,V2). 


Proof. Let 71 and 72 be two monotone paths that connect (ui,vi) with 
(r 4 ,(ui, ui), 0) and {u 2 ,V 2 ) with (r^(u 2 , U 2 ), 0) respectively. Let us assume that 
ri{ui,vi) > ri{u 2 ,V 2 )- As Figure |4] shows the paths 71 and 72 intersect at some 
point (uq, Vo). Therefore, a monotone path that connects (u 2 , 'C 2 ) with (r 4 ^(wi, tii), 0) 
exists. This path consists of a segment of 71 from (r^(ui, ui), 0) to (uo,vo) and a seg¬ 
ment of 72 from (uo,vo) to (^ 2 ,^ 2 ). This means that the point (u 2 ,V 2 ) is reachable 
from a point that is located to the right of the point (ri(u 2 , V 2 ), 0). This contradicts 
with the definition of function rj,. Therefore, the assumption ( 142 ,^ 2 ) 

is proved to be wrong and the first statement of the theorem is proved. The proof 
of the second statement is similar. □ 


Due to convexity of n Dij and monotonicity of functions and they 
satisfy the following recursive relations that will allow to use some sort of dynamic 
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programming for their computation. For {u,v) G TRij it holds that 


II 

(exjr] 

{uv') G g; n LB,j 

u < u,v' < J 

( 2 ) 

and for (u, v) G g^ Ci LBij 





-)• 

II 

extr 1 

{u'v') G g^ DTRij 

u' > u,v' > v'^j . 

(3) 


Relations m and ([3|) immediately result in the following lemma. 


Lemma 4. For each pair i G {1, ..., 2m}, j G {1, ..., n} the pointer r^ is constant 
on Pi n Rij and the pointer r^ is constant on fl Bij. 

Proof. For each pair of points (u, v) G Rij and (u', v') G LBij the condition u' <u 
in © is fulfilled and may be omitted. For each point {i,v) € 34,0 Rij a point 
(it', v') G n LBij exists such that v' < v. Consequently, this condition in ([2]) also 
may be omitted. Relation ([2]) becomes 

r 4 ,(u, v) = ri(^ extr {(u', v’) G giH LBij}^ , 

and ri{u,v) does not depend on u and v, which proves the first statement of the 
lemma. 

For each pair of points {u,v) G Bij and {u\v') G TRij the condition v' > v 
in Q is fulfilled and may be omitted. For each point {i,v) G 5 ^ H Bij a point 
{u',v') G g"^ DTRij exists such that u' > u. Consequently, this condition in ([3]) also 
may be omitted. Relation ([3|) becomes 

r^{u, v) = r"*'^extr {(m', v') G g^ Ci TRy , 

which proves the second statement of the lemma. □ 


Now relation m can be written in more detail. Let r* be the constant value of 
the pointer function r^ on g^ n Lij for some (i, j) and 





as it is shown on Fig[5a1 



If pj, n Bij = 0 then for (u, v) G 54 . C 



* 

II 


(4) 

If pi n Bij ^ 0 then for (m, j) G 5 ; fl and for 

{i,v) G 5 ; n Rij 


fr* 

ii u < a, 

(5) 

ri{uj)= < r;('U,j-l) 

a a < u < b, 

( 6 ) 

1 —1 

1 

ii u >b, 

(7) 

ri{i,v) = ri{b,j - 1 ). 


( 8 ) 
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(a) Pointer function on the upper bor¬ 
der depends on on the lower and left 
border of the cell. 


r* 



(b) Pointer function on the left border 
depends on on the upper and right bor¬ 
der of the cell. 


Fig. 5: Recursive dependency of pointer functions. 


Similarly, relation ([3]) can be specified. Let r* be the constant value of the pointer 
function on H Ty for some (t, j) and 



{i,v) G 


n R^J 


as it is shown on Fig l5bl 

If g'l' n Rij = 0 then for (u, v) G Ci LBij 




rt 

{u, v) = r*. 



(9) 

for {i 

- 1,-^) 

eg 

4 n Lij and 

for 

{u,j - 1) G gTl Bij 




1 

f r* 

if V 

>b, 

(10) 

- 

- 1,U) : 


r'^{i,v) 

if a 

VI 

V 

(11) 



1 

[ 

if V 

< a, 

(12) 

r^{u, 

j - 1) 

= r 

4(i,a). 



(13) 


The following Lemma [5] states that restriction of to g^ fl Rij is a piecewise 
constant function. Lemma [ 6 ] states that the restriction of ri to gi n Rj is also 
piecewise constant with the exception of at most one interval, where it is linear. 
For any set S we say that / is a partition of S' if S' = IJ int and int n int' = 0 

int^I 

for each pair int, int' G I, int ^ int'. 


Lemma 5. For each {i,j) a partition I^{i,j) of the set g^ H Rij to intervals ex¬ 
ists such that the function is constant on each interval in R(i,j): moreover, 
|/t(*,j)|<4m + l. 
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Proof. Indeed, either g'^ n R{ 2 m)j is empty or pointer is constant on 5 ^ n R( 2 m)j 
and equals 2m. Therefore, /^(2m, j) consists of no more than one interval. It follows 
from recursive relations dSD-dni) that there are no more than two intervals that 
belong to — l,j) and do not belong to The first interval comes from 

(uni). The second interval comes from (ini). The pointer r^{i — l,ti) is constant on 
each of these two intervals. Therefore, |/^(2 to, j)| < 1, \I^{i — 1, j)| < \I{i,j)\ + 2, 
and finally j)| < 4 to — 2i + \< 4to +1. □ 

Lemma 6. For each {i,j) a partition I\,{i^j) of H Ty into intervals exists with 
the following properties: 

- there is no more than one interval int G I{i,j) such that ri{u,j) = u on int; 

- function is constant on all other intervals; 

- moreover, \Ii(i,j)\ <2n + l. 

Proof. By definition, r 4 ^(u, 0) = u for each (u, 0) S n TiQ. Therefore, partition 
ll{i,0) consists of no more than one interval. It follows from recursive relations 
dl-dHl) that there are no more than two intervals that belong to Ii{i,j) and do 
not belong to Ii{i,j — 1). The first of them is included according to relation ([5]). 
The second interval appears in partition when u > 6 in relation ([7|). The 

function ri{u,j) is constant on each of these two intervals. Therefore, |/ 4 ,(z, 0 )| < 1, 
I4(b j)l < I4(b j -1)1+2, and finally \Iiii,j)\ < 2j + I < 2n + I. □ 

According to Lemma |4] the set g^ D B and the restriction of a function rl to 
this set can be expressed with subsets g^ D Bn, i S { 1 , 2 ,..., 2 m}, and values rj 
of a function on these subsets. According to Lemma IH the set gj, n T and the 
restriction of rj, to this set can be expressed with the sets I^ii, n) of intervals int and 
with numbers r™*, where int € Ii{i,n), i € {I, 2,..., 2m}. Numbers r™* determine 
the function on int so that if r}"* is less than the right endpoint of int then 
ri{u,n) = rff* for all (u,n) € int. Otherwise, ri(u,n) = u for all {u,n) G int. 

Lemma 7. Let X and Y be closed polygonal curves and for each i G {I, 2,..., m} 
the following data be known: 

- the set gl H Tio with pointer rj ; 

- partition I^{i + m, n) of g^ n 

- value for each interval int G I^ii + m,n); 
then testing S{X,Y) < e takes 0{mn) time. 

Proof. According to Lemmas [T] and [2l inequality 6{X,Y) < e is equivalent to the 
existence of a number u G [0, m] that fulfills ([T]). Such u exists if and only if a triple 
i G {1,2,..., m}, int G Ii{i + m, n), u G [0, m] exists that fulfills conditions 

(m, 0 ) e gl n Tio, {u + m,n)Gint, u + m<rj, m < r 4 ,(u + m, n). (14) 

The condition u < r^(u + m,n) in (jI4|l can be replaced with condition u < r}”* 
independently of whether the pointer rj,(u,n) takes constant value r}"* on int or 








May 18, 2015 0:39 WSPC/Guidelines HausdorfArticleEng2 


10 Schlesinger M.L, Vodolazskiy E.V., Yakovenko V.M. 

r|(M, n) = u. In both cases, condition (IT4l) is equivalent to 

(m, 0 ) G n TiO) {u rn,n) G int, u + to < rj, u < r™*, 
or in more detail 

Ci < u < di, — m < u < 5™* — to, u < rj — m, u < r]"*, (15) 

where and are the left and right endpoints of the interval int and Ci and di 
are the horizontal coordinates of the leftmost and the rightmost points of fl T^q. 
In turn, a triple i G {1,2,..., to}, int G Ii{i + to, n), u G [0,to] that satisfies (flSl) 
exists iff a pair i G {1,2,..., to}, int G Ii{i + rn, n) exists that satisfies 

max{ci, a™* — to} < min{di, 6 *"* — m,rj — to, r}"*}. 

Testing the last inequality takes constant time for any pair (i,int). According to 
Lemma [ 6 ] the number of intervals tested for each i does not exceed (2n + 1). There¬ 
fore, the number of tested pairs {i,int) is 0(mn). □ 

5. Computing pointer functions 
5.1. The general scheme 

The partition of T and B into intervals that represent and is obtained by two 
similar independent algorithms called the forward and backwards pass. Both passes 
consist of 2mn steps and in each step compute a pointer function on a border of 
some cell. From now on we will refer to the pair {i,j) as the number of the step. 

The forward pass proceeds cell by cell from left to right and from bottom to 
top and computes pointers on TRij based on pointers on LBij according to 
(l4])-([7l) starting from cell (1,1). The result of the forward pass are the partitions of 
T{i+m)ni 1 < i < TO, that represent on T. 

Similarly, the backward pass starts from cell {2m, n) and proceeds from right 
to left and from top to bottom and computes on LBij based on on TRij 
according to ©-(US]). The result of the backward pass are the constant values rj of 
the pointer function on n Rq for 1 < i < to. 

During the forward pass the function on gi n Tij as the partition I^{i,j) 
of Tij is stored in the following data structure. The function on each interval 
int G Ii{i, j) is represented by a triple end™*, r™‘), where 605 ™* and end"^^* 

are the endpoints of interval int and r™* has the meaning defined right before 
Lemma[3 The triples {beg, end, val) sorted by their endpoints are stored in a double- 
ended queue (deque) with the following operations performing in constant time: 

- reading and removing either the leftmost or the rightmost triple; 

- pushing a triple either to the left or to the right end of the deque. 

For a given number x the operations of cutting the deque to the left of x and 
cutting to the right of x are defined. Cutting to the left of x means removing 
all triples {beg, end, val) where end < x from the deque and replacing a triple 
{beg, end, val) where beg < x < end with a triple {x, end, val). Cutting to the right 
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of X means removing all triples {beg,end,val) where beg > x and replacing a triple 
{beg, end, val) where beg < x < end either with a triple {beg, x, val) if val < end, or 
with a triple {beg,x,x) if val = end. Since the triples in the deque are sorted, the 
time spent on cutting the deque is proportional to the number of triples removed 
from the deque. Note that each cut of the deque performs only one push to the 
deque. 

The partitions I^{i,j) of the sets g^CiRij and functions on g'^CiRij are stored 
in the same data structures. Of course, instead of the leftmost and the rightmost 
triple there are the lowest and the highest triple. 

The forward and the backward passes rely on the sets gi 0 Rij, gi 0 , g^ fl Rij 

and g^ 0 T^- to be precomputed. It has been mentioned at the beginning of Section 
m that these sets can be computed in a straightforward way. 


5.2. The forward pass 

The forward pass works with 2m deques Qi{i), i G {1, 2 ,..., 2m}, whose content 
depends on the number {i,j) of the step. The (z,})-th step starts with the known 
value r*{i — l,j) of on and with Q^{i) representing on Bij. The {i,j)-th 
step updates the deque Qi{i) to represent on and computes the constant 
value r*{i,j) of on Rij. 

As it has been shown during the proof of Lemma [ 6 ] a partition of non-empty set 
gi n Bii consists of a single interval. Therefore, at the start of the {i, l)-th step the 
deque Q\,{i) is either empty (if giCiBn = 0) or includes a single triple {beg, end, end) 
where beg and end are the horizontal coordinates of the leftmost and the rightmost 
points in giCiBn. The pointers r*{0,j) are obviously equal to 0 for each non-empty 
9l n Lij. 

When all four sets gi 0 Rij, g_i fl Rj, g^ fl Lij, g^ n Bij are non-empty the update 
of Qi{i) and the computation of r*{i,j) is done as follows. Let 


[a,6] = |M (u, j - 1) G 5^ nBy j , [c,d] = |M (w,j)G5tGry | 


1 . if c < a then push {c,a,r*{i — 1, j)) to the left end of Qi{i); 

2 . else cut Qi{i) to the left of c; 

3. read a triple {beg, end, r*) from the right and save the value r*; 

4. if 6 < d then push {b, d, r*) to the right end of Qi{i); 

5. else cut Qi{i) to right of d; 

6 . r*{i,j) = r*. 

The Operations 1, 4 and 6 represent relations ([5]), ([7]) and (|S]) directly. The relation 
(ini) is represented with Operations 2 and 4 indirectly in a sense that the uncut part 
of Q_i{i) remains unchanged. 

It is not necessary to consider all special cases when some of the sets gi 0 Rij, 
giDTij, g^DLij, g^DBij are empty. In each of these cases the {i,j)-th step consists 
of some part of Operations 1-6 or their slight modifications. Since we are mainly 
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interested in the computational complexity of the algorithm, considering the above 
described complete case is sufficient. 

One can see that one step of the forward pass does not perform in constant 
time due to the cutting of the deque that can take 0{j) time for the (i, j)-th step. 
Nevertheless, the complexity of the forward pass is 0{mn) as the following lemma 
states. 

Lemma 8. It takes 0{mn) time to complete the forward pass. 

Proof. The forward pass starts with initializing the deques Qi{i), i S 
{1,2,..., 2m}, and the pointers r*{0,j) = 0, j € {1,2,..., nj. Evidently, it takes 
0{m + n) time. 

No more than three triples are pushed in deque on each step. One triple is 
pushed either with Operation 1 or with Operation 2, the second is pushed either 
with Operation 4 or with Operation 5. The third push is made with Operation 3 
because the triple that was read and removed from the deque has to be returned to 
the deque. Therefore, no more than 6mn triples are pushed during all steps. 

Reading and removing triples are fulfilled with Operations 2, 3 and 4. Number 
of these operations in (i,j)-th step may differ from number of insertions in this 
step. However, the number of readings and removing during the whole forward pass 
cannot exceed the total number of insertions, consequently, cannot be greater than 
6 mn. □ 


5.3. The backward pass 


The backward pass works with n deques j € {l,2,...,n}, whose content 

depends on the number (*,j) of the step. The (i,j)-th step starts with the known 
value r*{i,j) of the pointer r^ on g'^ n Ty and with the deque Q^{j) representing 
r”! on n Rij. The (i, j)-th step updates Q^{j) to represent on g'^ n Lij and 
computes the value r*{i,j — 1 ) of the pointer r^ on g'^ n Bij. 

When all four sets g'^ fl Rij , g^ fl , g'^ n Lij , g^ n Bij are non-empty the update 
of Q^{j) is done as follows. Let 


[a, 6] = |z; {i,v) € g'^ n Rij'^ , [c,d] = ^v (i - 1, u) G g’’'fl Ly | 


1. if b < d then push {b,d,r*{i,j)) to the upper end of Q^{j); 

2 . else Q^{j) is cut to the up of d; 

3. read a triple {beg, end, r*) from the lower end and save the value r*; 

4. if c < a then push (c, a, r*) to the lower end of Q^O);; 

5. else cut Q^{j) down of c; 

6 . r*{i,j - 1) = r*. 

The Operations 1, 4 and 6 represent relations cni). dia and m directly. The 
relation (HU is represented with Operations 2 and 5 indirectly in a sense that the 
uncut part of Q^{j) remains unchanged. 


Lemma 9. It takes 0{mn) time to complete the backward pass. 
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Proof. The proof is based on the same considerations as the proof of Lemma [H □ 


6. The result 

Theorem 1. Let X and Y be closed polygonal curves with m and n vertices and 
S{X,Y) be the Frechet distance between them. Testing the inequality S{X,Y) < e 
takes 0{mn) time. 

Proof. Testing the inequality d(A, T) < e is reduced to the following computations. 
Computing the sets fl Ty and n Rij, 0 < i < 2m, 0 < j < n, is reduced to 
solving 2mn quadratic equations and takes 0{mn) time. 

Computing subsets g^ PI Tm and g^ n Bn, 1 < i < 2m, takes 0{mn) time. 

Due to Lemma [8] computing the restrictions of r^ to gi H takes 0{mn) time. 
Due to Lemma |9] computing the restrictions of r"^ to g^ H Bn takes 0{mn) time. 
Due to Lemma [7] testing the inequality 6{X,Y) < e based on these data takes 
0{mn) time. □ 
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